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VIBRATION AND CONTROL OF AXIALLY MOVING BELT SYSTEM 
(3 rd Report; Analysis by Parametria Excitation) 



Hiroki OKUBO, Kouetsu TAKANO, 
5 Osarai MATSUSHITA, Keiji WATANABE, 

and Yoshi HIRASE 



High speed operation is required for industrial machines 
equipped with belt driving systems. However, a belt has a 

10 limited speed due to occurrence of bending mode vibration 
resonance caused by machine errors such as eccentricity of 
pulleys. To overcome such resonance of a linear system, a 
parametric excitation method that provides fluctuation in 
belt tension is proposed through open- loop vibration control . 

15 This paper analyzes two types of fluctuation excitation: one 
is featured by frequency synchronized with pulley revolution , 
which is the most effective way, and the other is featured 
by frequency synchronized with two times the pulley revolution , 
which is an impractical way. According to analysis of the 

20 results of the well-tuned experiments, our control looks 
promising. 
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1* PREFACE 

Belt drives with rubber belts as power transmission 
means are widely used in various industrial units . Especially, 
running belt systems take a variety of forms according to 
5 the object and usage and are required various performances. 
High precision and high speed are required in running belts 
and resonance of lateral oscillation of the belts has thus 
become a major problem. As research with regard to oscillation 
restraining control, Chung et al. made theoretical analysis 

10 of active oscillation control by boundary displacement input . 
Also, Rahn et al. made theoretical analysis of oscillation 
control by parametric excitation. These are based on state 
feedback control and it is thus considered that there exist 
practical problems such as installation of sensors and types 

15 of actuators causing tension variation. 

The present research is directed to restraining 
lateral oscillation of a flexible means such as a rubber belt ♦ 
In the previous report, comparative study was made through 
basic experiments and analysis. In the present report, we 

20 are going to make a one -degree -of -freedom theoretical 
analysis of oscillation restraining of parametric excitation 
control and to clarify basic concept . We bear in mind a method 
that: is easy to practice, in which torque Ve±*.-.*.rtt;Um is given 
to a belt drive motor to cause tension fluctuation of a belt 

25 thereby restraining oscillation, as shown in Fig. 1. 
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Fig. 1 Driving Belt model 



General methods to control oscillation are 
10 absorbing oscillating energy through damping force to control 
oscillation, offsetting vibration waves by generating waves 
of inverted phase relative to the oscillation that has 
occurred, and the like . The former method is a way to feedback 
control the occurrence of oscillation. The latter method is 
15 away to open-loop control to excite so as to cancel oscillation 
that has occurred. Therefore, both methods differ in object 
and usage. In this report, we are going to make analysis of 
flexible or soft model for oscillation control of open-loop 
using parametric excitation. 

20 

2, SYSTEM MODEL 

As ii>eiri; Voaeti in the first revolt , a running bolt 
handled in this research has an extremely srr.all Ccricli'c 
25 force by the speed and has elgenf requency that t*?ill not vary 
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greatly according to the increase of speed. Therefore, we 
consider that influence of a speed term does not exist. Also, 
the source of excitation of the present belt system is an 
eccentricity of an idler pulley* According to the previous 
report, equation of motion in this case is as follows: 

i:4-a>«(l + Acos {vt + a)}x+2%a>ni 
— e*ti? cos cot+g* (1) 

wherein xi displacement, co : frequency, co „: 
resonance frequency, h: degree of modulation, v : excitation 
frequency, a : phase, t: time, e: eccentricity, g: 
gravitational acceleration, e*=(2/x)e, g*=(4/ k )g. C, e* , 
hi extremely small value, and co^co n . 

The second term on the left side designates tension 
variation and the third term damping. The first term on the 
right side designates eccentric centrifugal force. A number 
of reports were published regarding stability of March 
equation in the case of small non-linearity. 

The gist of this report is to consider whether 
resonance oscillation control is possible with the parameters 
of {h) and ( a) for tension variation. Non-linear osoi llation 
analysis ic made using asymptote expansion method . In urder 
to make a later analysis easy, the differential equation is 
rewritten as follows : 
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X -f (0 2 X = (CD 2 ~ Q)1)X ~ 2 %G)n± H ' CO 2 COS COt 

-G> 2 nkcos(vi-\-a)z+g* (2) 



3. OSCILLATION ANALYSIS FOR FLEXIBLE SYSTEM 

3.1 Arrangement of Differential Equation 

The following gives an analogous solution to 
equation (2) . 

x=a cos (o)t + 9) + eui(a } (9, ut) (3) 

wherein a and 0 designate amplitude and phase, 
respectively, which are not constants but variables that vary 
slowly due to non-linearity. They are supposed to be as 
follows : 



^=eAx(a, d) + ^Ma> d) 



Trie left side of equation ( 2 ) ic rewritten as follows s 



jf + co 2 x — e(^—2A\w sin $ 

[{A^^B x ^-aBl-2acoB^ cos <l> 



+ e 2i 
+ 



1 dd 



) j sin ^ 



+ 2A 



dadt^ 1 * 1 dddt] 



10 



(4) 



In this paper, we consider the case of v=co where 
excitation is carried out at the same frequency as resonance 
frequency and the case of v=2 co where excitation is carried 
15 out at the doubled frequency. When ignoring the small second 
order term and arranging the right side of equation (2) 
to obtain the following: In this case, <p = cot+ 6 . 



(<y z -<y«)x=(<y 2 — Q>l)a cos 0 
-2K<D n ± =2 %o)nam sin <f> 

+e\ Ai cos <p-aBi sin 

e*«f cos cut = e*of{co5 d cos ^-rsm </ sin 4>) 
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At v = co g 



— colli cos ( vt 4- a)ar 

-hcos(a-(9)} 

— a) 2 heu\ ccs + a) 



At v =2co , 

10 

— ait/i cos (y/ + ff):c 

---^^cos(a-2<9) cos<* 

+cos (a-2$) cos 3^ 
-sin (a-2d) sin 3^ 
15 — sin(ff— 2(J)sin 0} 



Both sides of equation (2) are made equal to each 
other using the above-mentioned equations. 

20 

3.2 Analysis of Amplitude and Phase (at v =2co ) 

In the case where the excitation frequency is twice 
as hiqh as the frequency, or v**2to , by tiOmparxng cos with 
sin ^ in the equation that was acquired by ignoring the terms 
25 more than a 2 , the following equations can be obtained: 
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da y e*w . 
—~- {wM ~ -sin 



Am 



a sin (a-2d) 



dd _ to 2 - col e*co ^ o 
-df- 2^ 2T ccsd 

+^cos(a-2d) 



(5) 



10 uj will be obtained from direct current and harmonic 

component that are not coefficient in the basic wave of cos 
</> and sin (/> - 



/^4-,^ W— ^{cos ( a -2d) cos 30 

15 

-sin(a-2<?)sin3^} (6) 



Therefore, e u 2 will be given as a special solution 
to the above equation. 
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eu\ =-~r+ Ci cos 3^+ C 2 sin 30 ( 7 ) 

Cz= —fz-sin'sCt-zo) 
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Here, new variable u(t)= a cos 6 , v(t)~ -a sin d 
are introduced. Then, a and 0 are rewritten into the 
following: 



a 2 =u 2 +v 2 



(9 = arc -tan — 

u. 



•(8) 



Using these equations, equation (5) will be 
10 rewritten as follows: 



[ft Q* 



-(9) 



Wherein, 



15 



20 



02 = 



co 2 — col , colh 



2<o 



H — r^cos a 
4co 



# 2= ~ co 2 -col j o)\h 



1® 



4(0 



cos a 



J. J /.^aiysis af Amplitude and Phase (at v^oj) 

In the case where the excitation frequency equals 
to the f rpqurnny . or v = t by comparing coefficients of basic 
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waves of cos 0 and sin $ with each other in the equation 
that was acquired by ignoring the terms more than f 2 , the 
following equations can be obtained; 



da 

di ' 



• — Zcona 2~sm d 



dd _ ay 
dt ~ 



2 

2^ U~ 003 * 



•(10) 



£ Uj will be obtained as a special solution to 
10 non-basic waves. 



+ cos(a-(?)} (11) 



15 

Therefore, £ u, will be given as follows: 



20 



cos (20 + a- <9) • 



(12) 



When the excitation frequency is equal to the 

frequency, a tcrir. of s:iodt2lat;lori degree h i n ?soi: included in 
the analysis up to a minute term of degree one, as can be 
25 seen in equation (10). Therefore, effect of 7 parametric 
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excitation cannot be found. Consequently, analysis including 
minute terms of up to degree two will be carried cut 0 The 
minute term of degree two at the first term on the right side 
of equation (2) will be as follows; 

5 

• e(a> 2 - dQm =W - cos <* - *> 

+^ccs(2^+s-<9)} 

The minute term of degree two at the second term 
10 on the right side of equation (2) will be as follows: 

-2^(-fcos^-^sin^^) 

The minute term of degree two does not exist at the 
15 third term on the right side of equation (2) . The minute terra 
of degree two at the fourth term on the right side of equation 
(2) will be as follows: 

- wlh cos (cot + a)e U\ 

Xcos (a— d) cos <l> 

+ (*)'*k-Hr^(«-*>} 

Xcun (a 3) sin ? 

-(it)"*'* 03 * ■co3(3 5 >+2ar-2S>} 

25 
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da/dt, d0/dt will be obtained as follows; 



10 



-sin 8 



da 



d& ^ a> 2 -col (a? -ml) 2 
dt 2<o 8<y 3 



2a> 3 



sin 



2a 



-cos (9 



(13) 



Using u(t) andv(t), equation (13) is rewritten into 
the following: 



15 



20 



j:W2 ?t:w;i <»> 



Wherein , 



/ 5 . = -^„+-^-sin2a 

^' 8a» s + "l25 r ~ "W 008 2ff 

R, ~ 2^ 8a, 8 

tt/ — Wnft-y | T/ Cutirig 



r~sxz jy 



25 
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4= OSCILLATION ANALYSIS AT RESONANCE 



In this paper , we are going to proceed with 



oscillation analysis at the time of resonance. At resonance, 
co -co n is satisfied . Then, we are going to obtain resonance 
5 oscillation at steady state in the case of v = co and v =2 
co of parametric excitation and discuss stability of the 
steady state. Finally, we seek the condition in which 
reduction in oscillation is stably possible by parametric 
oscillation . 



greater than the frequency at resonance can be obtained by 
substituting ^for co , 0 for du/dt. and 0 for dv/dt in equation 
15 (9). Then, P 2n , Q 2n , R 2n , S 2n are substituted for B 2 , Q 2 , R 2 , 
S 2 , respectively, u and v are shown as follows: 



10 



4.1 Amplitude at Steady State (at v-2co) 



Steady state where excitation frequency is twice 



Qzn e*0)n 



V 




(15) 



20 



As the amplitude 1s shown by a 




25 



is obtained. Without control, h equals to 0, 
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Amplitude a p0 at resonance without control is shown by 



5 That stands for amplitude for resonance oscillation 

in linear system, which is obtained by multiplying 
eccentricity by resonance magnification, 

4.2 Reduction in Oscillation and Stability at v=2<o 
10 At the time of resonance, co n is substituted for co 

in equation (9) . There are two case of h=4 £ and h=£4t> . Also, 
initial state of u and v is set at u(0)=u o , v(0) = vb/ 
respectively. 

15 9At h~4 £ 



20 



(18) 



vit^Av + Bvt + Cve- 2 ***'* 
wherein , 



A u = g *g^ a +y(ao(l-sin a) + t*cos a) 
i 

5 U ~ Ct};z COS Ot 



H — =-t?i*> cos a + t/oU+sm a;; 



14 



B v = — r£*Wl+sin or) 



Co=— -sin a) 
— |-(ho cqs a - vo{\ -sin or)} 



10 



In order that amplitude a does not diverge , equation , 
15 B u =BvX=0 should be satisfied. At this time, a should satisfy 
the following condition: 

f (19) 

20 In this case , a does not converge if a is even slightly 

out of the condition of equation (19). However, in actual 
machines , it is difficult to set a without any error . Therefore , 

fr~4 C is not a practical case. 
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10 



zt(t)=Au+B tl e-°"«+'" t " + C u e- 61 * l< -'" iU 



(20) 



wherein , 



a _ 2g* fecos a 
^ u ~~ 2,2 _i a** 



/H-iT" — 2"U-t-sin a) — g~ cos a 

Ctt =-^^^(i-sina)+-f cos a 

A - 2g^4 f+Asing) 
A " V-16? 2 

o _ e*(l— sin a) wo ^ _ , t>o/i • \ 
°" fe+4g — 2~ 2^ ^ 



15 



C 0 = h—A% ~' "i — cos ff-h^Ufsm 



20 

In order that amplitude a does not diverge, h<4C 
should be satisfied so that the third term of equation (20) 

can dduit). 

Next, Fig. 2 shows amplitude a at the time of 
2* convergence. Horizontal axis of Fig, 2 designates phase angle 
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a . In addition, amplitude at h^4 C is the one when the initial 
value is set at u 0 -0 and v 0 =0. Also, in Fig. 2, vertical axis 
a=e*/2 £ designates resonance amplitude at the time of 
linearisation. In this parametric excitation control, when 
5 a<e*/2 C can be achieved it contributes the reduction in 
oscillation. It can be found that when a~-n/2 is selected 
efficiency becomes better. Also, even in the case where best 
selection was made, a is up to a^e*f 4 C • It can be seen that 
reduction in oscillation is limited to a half of the 
10 oscillation. 



15 



20 
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Fig* 7, Amplitude at vibration resonance by parametric 
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4.3 Amplitude at Steady State at v^co 

Steady state where excitation frequency at 
resonance Is equal to the frequency can be obtained by 
substituting co n for & , 0 for du/dt, and 0 for dv/dt in equation 
5 (14) . Then, P ln , Q ln , R lD , s ln , W UD , and W VD are substituted for 
Pi* Qi, Ri, Sj, W u , and w Vt respectively, u and v are shown 
as follows ; 



10 



SlnWun-QlnWvn 
P\nS\n~ Q\nR\n 

~RluWun + F\nWvn 
PlnS%n~~ QlnRln 



(21) 



Amplitude a p can be calculated by equation ( 8 ) . 
Amplitude a p0 without control is similar to equation (17), 

15 

4.4 Reduction in Oscillation and Stability at v = co 

At the time of resonance, co n is substituted for co 
in equation ( 14) . Then, u(t) and v(t) will be obtained. Also, 
initial state of uand vis setat u(0) = u 0 , v(0) = v 0 , respectively. 
20 There are two case of h 2 =24t,/4s and h 2 *24t> / 45 . 



mat h*°24{/& 



25 A ""2ft»n '4f4s- Ba 2£m„ 
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15 



20 



Wherein , 



Du — Uo 

En— K<i)nitti{\ +~^ r sin 2aj 

+-^=«o(2- 3 cos 2a) Hp^- sin a 
F„= -ipU 3 V sin a-^-n* cos <r 

-^<2-3cos2a) 



Ey . Fy 



Wherein, 



?5 
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E 0 = — ^*«,(2+3 cos 2ff) 

+ %C0 n v o (l — sin 2a) 
■ -/gfrg* g*a>« 

+ -5 T7 vr C0SC — ~ 

f,=^rr/ cos a--jp^V sin a 



(l--^ sin 2a) 



In order that amplitude a does not diverge, B u = fl v =0 
should be satisfied. Then, A will be as follows: 

10 

j e*a>l 2-3 cos 2a (23) 

g* /5"sin or + cos a 



•At h 2 ^24^1 41 

15 



u(t)=A u +Bue- p, + C a e-'" 
v(t)=Aa+B B e- pt + C B e-" t 



20 



} (24) 



wnsrein. 



25 
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10 



~ <, , ojJi 2 Vz hg* . 

+ Wr g" < sin 2ff-t/dcos 2a) 
F a =--^f cos a--^ sin a 
-^ffV 3 cos 2a) 



Wherein, 



£) 0 =fo 



^■iuo cos 2« + w> sin 2a) 



8 

2 + 2a>„ cos a 
F„=-^cosa-^f sin a 

20 +— jg—sin2ff 2^ 



[n order tor amplitude a not. to diverge, 
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should be satisfied. 

Here, a new parameter is introduced . ( can be 
rewritten into the following using e*. co n7 g* g and Z- 



?=z(-^#-) 2 (26) 



Z takes a positive value, which is specific to a 
belt. As the value becomes greater, the better damping 
properties the belt has. Next, we are going to find h that 
10 makes amplitude zero. Here, by substituting KC for h 2 , 
analysis will be carried out. That is, 

wherein, K is a positive constant. 
15 On the other hand, in order for amplitude to have 

a convergent value, from the condition of equation (25), 



«<7r (28) 



20 F a 2 + F V 2 =Q is required to achieve a=0. F u 2 +F v 2 will 

be rewritten into the following: 

Fi+I<i=FcZ+fi l tZ ir ~+Fc (29) 

wherein , 

25 
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F a - 



K* , K 
23C4 ' 4 



_ IK" 2 



K m 



ccs a— "g sin g+ "2 cos 
F c =-~ ^-(13-12 cos 2ff)-y|sin 2a-i~j 



') 



Equation (29) can be considered a quadratic 



expression of z 



1/2 



In order to make amplitude zero, the 



following condition is required: 



10 



4F< 



•(30) 



15 At this time , the condition of a and K that satisfies 

equation (30) is as follows: 

* = ° } (31) 

KZ=V 

20 From equation (28) , the following is a belt condition 

that can make amplitude zero: 

y->/ 5 r-n t\w . (32) 

25 
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10 



15 



Here . we assume a belt system of Z-l • 28 . Fig . 3 shows 

an amplitude at resonance oscillation In the case of h 2 = £ , 
h 2 =2 C , h 2 '~*4 Z, . That is, in this case, it can be seen that a=Q 
is possible at <z=0 and h p±T : 0 * 78. 



2.5 
2 

1.5 
1 

0.5 



2=1.28 



» 4 0 it 



o h l =0.78t 



-0.5 



0 
a 



0.5 



1 

xn 



Fig. 3 Amplitude at vibration resonance by parametric 
excitation (u=co~a) a ). Z=1.28 

5. OSCILLATION SIMULATION 
Numerical calculation was carried out in equation 
of motion (2) in order to search variation of x relative to 
h and a. Parameters used here are shown in Table 1. Initial 
state is x=e*/ C and x=0.0 m/s . 



20 



Table 1 Meaning of parameters 



-1- 



j 0,05 j 2$[Hz} I 25jff *j j 25;//<r) \^.l[mrn\ j 1,28 j 



25 
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10 



15 



20 



Fig. 4 shows envelopes of oscillation amplitudes 
of belt oscillation k relative to time in the case of v 
~2co and v-co. 

• At v~2cd, a=0 cannot be obtained. a^e*/4£ is the limit, 

• Since a belt with a higher damping properties of £=1.28 
is used, at v = co , it is found that a=0 can be obtained 
with adjustment of h and a. Also, it converges in about 
Is. 




(1) v=2a>h=3? o=-Jt/2 

(2 ) v=2coh=3.8C a=-n/2 



, \<2) <*) v=C0 =0.78^ 0t=0 



Oil- 



§5 



012345678 
t [sj 

Fig. 4 Amplitude Cuves of belt vibration (x) 
6. CONCLUSION 

When we consider reduction control of resonance 
oscillation of a belt via parametric excitation, we can 
conclude as follows? 

( 1 ) when excitation fregpuBnay Is twice greater than 
the frequency or v~2o>\ 
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In fact, since errors of a could be introduced, if 
we take a value of h that is smaller than 4 £ and that is 
as close to 4 £ as possible around a = ~K/2, reduction in 
oscillation can be achieved up to approximately half of it. 
5 (2) When excitation frequency Is equal to the 

frequency or i/=CO; 

Generally, it has been found that oscillation can 
be reduced when excited at tf=0 and in the stable range of 
• Especially, existence of belt condition (Z) was 
10 shown that can make amplitude zero at the time of resonance. 
It was also shown that excitation condition of perfect 
cancellation of oscillation is tf=0 and K-l/Z. 
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